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Abstract 

The effect of the Gauss-Bonnet term on the SU(2) non-Abelian regular stringy 
sphaleron solutions is studied within the non-perturbative treatment. It is found that 
the existence of regular solutions depends crucially on the value of the numerical factor 
(3 in front of the Gauss-Bonnet term in the four-dimensional effective action. Numerical 
solutions are constructed in the N = 1,2,3 cases for different (3 below certain critical 
values (3n which decrease with growing N (N being the number of nodes of the Yang- 
Mills function) . It is proved that for any static spherically symmetric asymptotically flat 
regular solution the ADM mass is exactly equal to the dilaton charge. No solutions were 
found for j3 above critical values, in particular, for f3 = 1. 



Since the Bartnik and McKinnon's discovery |l[ of the regular particle-like solutions 
to the coupled system of the Einstein- Yang-Mills (EYM) equations there was a growing 
interest in revealing their possible physical significance. It was shown that these solu- 
tion could play at the ultramicroscopic distances a role analogous to that of electroweak 
sphalerons @. Sphaleron interpretation is supported by the existence of the odd-parity 
YM negative modes fU (apart from the previously known even-parity ones 0]), as well as 
the fermion zero modes and the level-crossing phenomenon |J. Natural question arises 
whether the EYM sphalerons survive in more sophysticated field models suggested by 
the theory of superstrings. It was shown recently that regular sphaleron solutions exist 
within the context of the Einstein- Yang-Mills-Dilaton (EYMD) theory @, [0, f|, PL 



0. Remarkably, they have a dilaton charge exactly equal to the ADM mass. This prop- 
erty is similar to that of the extremal dilaton black holes which are likely (at least some) 
to represent exact solutions of the string theory. 

To further investigate possible relevance of the EYMD sphalerons to the string theory 
we study here the EYMD system with the Gauss-Bonnet (GB) term which is typically 
present in stringy gravity as the lowest order curvature correction. Similar problem for 



Abelian dilatonic black holes was studied recently within the perturbative approach fll 
However, for regular solutions, one needs a more precise treatment. In order to see in a 
continuous way how EYMD solutions are modified by the GB term we introduce into the 
lagrangian a numerical factor (3 so that — corresponds to the pure EYMD system. 
It turns out that series expansion of the regular solution near the origin is essentially 
/3-dependent. Also, computing the GB contribution into the effective energy density 
on the background EYMD solutions, one can observe that the GB effect becomes non- 
small for f3 of the order of unity. For this reason we avoid any perturbative treatment 
of the GB term and attack the problem numerically. Starting with (3 = we increase 
gradually the value of this parameter and search (using the shooting strategy) for solutions 
interpolating smoothly between the regular asymptotic expansion near the origin and 
an asymptotically flat expansion at infinity. Although the leading terms of expansions 
near infinity are not modified by the GB corrections, those near the origin are affected 
substantially. We construct numerical solutions for N = 1,2,3 and some (3 ^ and 
show that regular solution cease to exist above certain critical values (3^ depending of 
the number of nodes iV of the YM function. For all solutions found within the domains 
of existence, modifications due to GB term are relatively small, and all characteristic 
functions still preserve the typical behaviour they have in the pure EYMD case. We also 
prove analytically that the dilaton charge of any regular solution (with an exact account 
for the GB term) is equal to its ADM mass independently on the value of (3. For (3 = a 
stronger relation holds between g 0Q and the dilaton factor everywhere. 

We start with the following bosonic part of the heterotic string effective action in four 
dimensions in the Einstein frame : 

S = tL I {(- R + 2d ^ d ^)- a ' e M-2®)(F a ^Fr-f3G)}^d 4 x, (1) 

where <3> is the dilaton, F is the Yang-Mills field strength and G is the Gauss-Bonnet 
term which can be presented as the divergence of the topological current 

g = r^xtR^ - tR^nr + R 2 = V M tf" . (2) 
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Integrating by parts the GB term in (1) one can rewrite the action in somewhat simpler 
form (both in (1) and (2) we ignore surface terms which are not relevant for the present 
analysis) 

S = IcL / { + 2d ^ d "^ - a ' e ~ 2 *( F ^ F a V ~ 2/3(<yt>)i^)} V=9d A x . (3) 
We parametrize the metric of the static spherically symmetric spacetime as 

ds 2 = Wdt 2 - ^- - R 2 (d6 2 + sin 2 9d<J) 2 ) , (4) 

where W = wa 2 and all functions depend on the single variable r. In this case only the 
radial component of the topological current is relevant 

4(wa 2 )'(wR' 2 - 1) 
K = RW ' (5) 

Here and below primes mean derivatives with respect to r. 

A magnetic part of the static spherically symmetric SU(2) Yang-Mills connection can 
be expressed in terms of the single function of the radial variable /(r) 

AlT a dx» = (/ - 1)( VZ0 - L e sin 9d4) , (6) 

where L r = T a n a , L e = dgL r , = (sin 6')™ 1 (9^L r , n a = (sin 9 cos (f>, sin 6 sin <f>, cos 6) is 
the unit vector and T a are normalized Hermitean generators of the SU(2) group. 

Integrating out the angular variables in (2) and eliminating some total derivatives one 
obtains the following reduced effective action 

S = J dtdr (L g + L m + L GB ) , (7) 



where 



is the gravitational part, 



is the matter part, 



L g = | (R'(wR)' + 1) + wa'RR' , (8) 



L m = -\wRV - ^Fe^ , (9) 



L GB = 2a'po- l &W\wR' 2 - l) e - 2 * , (10) 
is the Gauss-Bonnet contribution, and 

(1 - f 2 ) 2 

F = 2wf 2 + [ -^^. (11) 

Note that an arbitrary rescaling of the slope parameter a' — > ka' together with the 
corresponding rescaling of the radial variable r — > \f~kr is a symmetry transformation 
of the effective action (7). Choosing Planck units a' = 1 we are left with the only 
dimensionless parameter (3. 
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The equations of motion (including an Einstein constraint) can be obtained by direct 
variation of (7) over a, w, R, /, $. Then fixing the gravitational gauge as R = r one finds 
the following set of equations 

•L = + 2 J^1 + - _ , (12) 

cr r r v a a 2 ' 

u> 1 H e + rw<P = (1 — 4pw(e ) J , (13) 




+ „ L r ^ _ ii^iM + 4/5 (^'^V = , (14) 



(w/e -)' + fZ(WV^ = 0, (15) 

(<7rW)' + ff Fe"» + 2,3 ^'(^-^ V^ = „ (16) 

It is useful also to compute an effective energy density as it enters the standard Einstein 
equations with account for the GB term 



2T ° = w& 2 + -^e~ 2 * + % {2w(w - !)($" - 2$' 2 ) + $V(3w - 1)} e" 2 * . (17) 



For (3 = the system reduces to that of @ and the corresponding solutions exibit typical 
BK structure of the YM function: solutions start from / = ±1 and goes asymptotically 
to =1=1 either monotonically (N = 1) or after N — 1 oscillations around zero. 

As a first step of the analysis we calculate the GB term and the corresponding density 
y/—gG substituting the sphaleron solutions found without an account for the GB term. 
Numerical results are shown on the Fig. 1. One can see that the value of GB term 
increases with growing number of nodes of the YM function. It can be anticipated that its 
influence on the sphaleron solutions will increase for higher N. We have also calculated 
the effective energy density (17) for the background EYMD solutions. Fig. 2 clearly 
shows that relative contribution of the GB term for [3 = 1 is not small. This presumably 
invalidate any attempt to treat the GB term perturbatively, so we are faced with the 
problem of constructing numerical solutions to the system (12)-(16). 

To define the ADM mass M and the dilaton charge D one writes asymptotic expansions 
for W 

2M 2D 2 M 

W = l ~- ~ + °(^)> ( 18 ) 

and the dilaton 

D DM 8M 2 D-D 3 1, 
® = ^ + J + — + ^ + °(^)- ( 19 ) 
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The corresponding expansion of a reads 

D 2 AD 2 M 

"=i-S5— ST + O^). (20) 

To ensure asymptotic flatness it is sufficient (as it is for (3 — 0) to have for the Yang- 
Mills function / 

f = ±l + 0(-). (21) 
r 

Clearly, GB-induced terms do not influence the leading behaviour of solutions near infin- 
ity. 

In contrary, an expansion of regular solutions near the origin is affected by the curva- 
ture terms. From the system (12)-(16) one finds 

/ = _i + fe r 2 + (r 4 ), (22) 

$ = $ + $ 2 r 2 + 0(r 4 ) , (23) 

a = a + a 2 r 2 + 0(r 4 ) , (24) 

W = W + W 2 r 2 + 0(r 4 ) , (25) 

or in terms of w: 

w = 1 + w 2 r 2 + 0(r 4 ) , (26) 
where the following relations hold 

W = 0-q, W 2 = 2a a 2 + w 2 a 2 . (27) 

Let us prove that for any regular solution to the system (12)-(16) (if exists), the ADM 
mass M is exactly equal to the dilaton charge D. Combining Eqs. (12), (14) and (16), 
after some rearrangment one can find the following identity 

( W'r 2 \ ' 

[2crrW+- =W, (28) 

where 

Q = a' 1 {(w - 1)(W' + 2W$') - 2r$wW} . (29) 
Integrating this relation over the semiaxis with account for (18)-(20) on gets 



W W'r 2 

— r 2 $' H 

a 2a 



= M - D = 2(3 [Q(oo) - Q(0)] . (30) 



Now from the expansions (18)— (21) and (22)-(26) it can be found that both above bound- 
ary values of Q are equal to zero, what proves the exact equality M = D. Remarkably, 
this property of regular EYMD solutions observed first in 0, remains true with account 
for the GB term for any value of (3. There is an important difference, however. In the 
case (3 — a stronger identity 

W = exp(-2$) (31) 
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holds, which is similar to the well-known relation for the extremal magnetic dilatonic 
black holes, where it ensures regularity of the metric in the string frame. When GB term 
is taken into account this is no longer true while the relation M = D exibiting the validity 
of (31) in the asymptotic region still holds. 

Similarly to the system of Einstein- Yang-Mills-Dilaton equations |j| , [[| , [§] , [[| with- 
out Gauss-Bonnet term there are three independent parameters in the series solutions of 
the system (12-16) near the origin: b, $o an d cr . From them the quantity $o is some- 
what trivial because of the symmetry of the system under a dilaton shift accompanied by 
suitable rescaling of the radial coordinate (if desired, exp(— 2$o) ma y be absorbed into 
redefinition of parameters in (22)-(26)). However, there is a substantial complication as 
compared with the pure EYMD theory. In order to fulfil the system (12-16) in the first 
leading order, the coefficient $2 has to be one of the real roots of the following algebraic 
equation of the forth order 

($ 2 + 2&V 2 * ) (l + 16/M> 2 e- 2 * ) 3 + 32/% 4 e- 6 * (l + 8/3$ 2 e- 2 * ) = 0. (32) 
Once $2 is found, two other coefficients W2 and 01 can be obtained as 

W2 = -l + 16^ 2 e-»o> ^ 



q e- 2# °(4fe 2 + 4/ 3w 2 $ 2 ) 
1 + 16/3$ 2 e- 



^ = ; , 7^ - 2 g n • ( 34 ) 



It is convenient to regard the Eq. (32) as giving the value of $2 a s a function of b, 
while $0 is fixed. In fact, a dilaton shift 

$0 -> $0 + (35) 

leads to a solution related with the initial one by a radial rescaling. Physically the 
normalization $oo = is preferable since it ensures a unique mass scale for all solutions. 
But technically is is convenient to solve the system first by fixing $0 arbitrarily, say, 
$o = 0. Then the rescaled solution will result from 

b -> 6exp(25$ ), $2 -> $ 2 exp(25$ ), <r -> °o- (36) 

At the final stage of the calculation we rescaled solutions imposing the condition = 
in order to fix a unique mass scale for all of them. 

The numerical strategy consists in solving the system (12)-(16) starting from the 
series solution (22)-(26) near the origin. The crucial role is played by the parameter b 
which shoul take a discrete sequence of values. For (3 = the solution of Eq. (32) reads 
$2 — — 26 2 exp(— 2$ ), an d clearly this does not impose any restriction on this parameter. 
But for (3 7^ it turns out that real solutions for $ 2 do not exist in some region of b. 
Hence, in addition to the problem of "quantization" of b one has to ensure that b belongs 
to region where the real roots of the Eq. (32) exist. It happens that if (3 is greater than 
some (iV-dependent) critical value @n, the allowed region of b does not contain those 
quantized values for which regular solutions exist. Only for /3 < (3 N regular solutions 
exist and exibit behaviour similar to that of the EYMD solutions. 
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Real roots of the algebraic equation (32) form two branches as shown on the Fig. 
3a,b in terms of the quantities b = 6exp(— 2$ ), ^2 = $2 exp(— 2$ )- F° r roots from the 
second branch (Fig. 3b) we didn't find any solution, they seem to correspond to b outside 
the above quantization domain. Note that this branch does not contain the EYMD root 
corresponding to (3 = 0. The first branch 3a has a solution for (3 = 0, while for any (3^0 
there are two negative solutions with absolute values $™ aa: (&) and $™ m (&). From these 
two, it is just the second one, <3>™ m (6), which has the limiting value $ 2 = —2b 2 exp(— 2$ ) 
when (3 — > 0. No regular solutions to the system (12)-(16) corresponding to $ max (fr) were 
found neither. 

Starting with the known ft = 0, N = 1 EYMD solution we increased gradually the 
value of (3 searching for the desired quantized b related to $™ m (fe). Numerical integration 
of the system (12)-(16) was done using the Runge-Kutta fourth order scheme. The 
values of the parameters for N = 1 case, found numerically for some (3 together with the 
corresponding ADM mass M = D are given in the Table 1. The solutions were rescaled 
to ensure $00 = 0. 



Table 1. N=l. 



p 


b 


$0 


O"0 






M = D 


w 2 


0. 


1.073 


0.9311 


0.3936 


-0.3576 




0.578 


-0.7153 


0.1000 


1.026 


0.9199 


0.3840 


-0.3523 


-3.390 


0.573 


-0.7344 


0.2000 


0.9866 


0.9122 


0.3744 


-0.3566 


-1.475 


0.568 


-0.7697 


0.3000 


0.9619 


0.9120 


0.3597 


-0.3833 


-0.8376 


0.563 


-0.8496 


0.3700 


0.9657 


0.9231 


0.3421 


-0.4938 


-0.5198 


0.560 


-1.0933 



One can observe that with increasing (3 two real roots $™ n (6) and <&™ aa: (&) converge 
and merge together for a limiting value (3\ approximately equal to 0.37. For (3 > (3\ there 
are no such b which could generate asymptotically flat solutions with N = 1 compatible 
with the existence of the real root $2(&) of the Eq. 32. 

Similar situation was encountered for higher-iV solutions. Numerical results for N = 2 
and N = 3 are presented in the Tables 2, 3. Figures 4-8 depict the corresponding 
numerical curves for some values of (3 and N. 



Table 2. N=2. 



(3 


b 


$0 






Qmax 


M = D 


w 2 


0. 


8.3612 


1.7923 


0.1665 


-3.8796 




0.685 


-7.760 


0.1000 


7.1902 


1.7481 


0.1529 


-3.5165 


-15.982 


0.673 


-7.558 


0.2000 


6.4017 


1.7297 


0.1370 


-3.6597 


-5.7461 


0.660 


-8.161 


0.2208 


6.3344 


1.7343 


0.1320 


-4.2904 


-4.3127 


0.657 


-9.478 



Note, that the numerical values of ADM mass/dilaton charge monotonically decrease 
with growing j3 for each fixed N. Also, it can be observed that the limiting values of 
(3n decrease with the increasing iV: (j3x = 0.37, (3% = 0.22, /3 3 = 0.21,...). It can be 
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Table 3. N=3. 



p 


b 


$ 


O"0 




Qrnax 


M = D 


w 2 


0. 

0.2000 
0.21117 


53.8351 
36.5453 
36.2683 


2.6920 
2.5930 
2.5956 


0.0678 
0.0504 
0.0492 


-26.600 
-22.348 
-25.151 


-30.543 
-25.212 


0.7042 
0.6744 
0.6726 


-53.202 
-49.817 
-55.558 



anticipated that (3n has a limiting value (3^ as A" — > oo, which would presumably give 
an absolut bound of the existence of static spherically symmetric regular EYMD-Gauss- 
Bonnet solutions. It is also interesting to note, that although the contribution of the GB 
terms to the energy density for (3 close to (3 N is not small (as it is shown on Fig. 5), the 
behaviour of / and metric functions is very similar to that of pure EYMD solutions. It 
has also to be noted that, when the limiting value (3n is approached, neither singularities 
no other numerical problems arise; so the only reason for the absence solutions when (3 
exceeds the above critaical value is an intrinsic incompatibility of the series expansion 
near the origin. 

We conclude with the following remarks. When Gauss-Bonnet term is included, the 
total number of derivatives in the system of equations increases, as well as the dergee of 
its non-linearity. However, in a limited region of the numerical factor (3 the behaviour 
of solutions remains qualitatively the same as in the pure EYMD case. Moreover, the 
remarkable equality of the ADM mass to the dilaton charge remains unaffected by the 
GB term for any (3. However it is likely that EYMD sphalerons are destroyed by the 
Gauss-Bonnet term for sufficiently large values of (3. The most persistent is the N — 1 
solution, which exists up to (3 = 0.37. Higher A" solutions cease to exist for lower f3, the 
limiting value is likely to be of the order of 0.2. 

This work was supported in part by the ISF Grant M79000 and by the Russian Foun- 
dation for Fundamental Research Grant 93-02-16977. 
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Figure Captions 



Fig. 1. GB term (B) and GB density (A), calculated for pure EYMD N = 1 solutions 
@, curves (C) and (D) - GB density for N = 2,3 EYMD solutions. 

Fig. 2. Contributions to the energy density r 2 * Tq , from YMD (f3 = 0) and GB parts 
{(3 = 1) calculated using EYMD solutions : (A): N = 1, YMD; (B): N = 1, GB; (C): 
N = 2, YMD; (D): N = 2, GB. 

Fig. 3a,b. Real roots of Eq. 32 (two different branches), (3 = 0.1,0.2,0.37,0.5,1 in 
terms of b = 6exp(-2$ ), $ 2 = $ 2 exp(-2$ ). 

Fig. 4. "Gauss-Bonnet" mass distribution (contribution to ADM mass from (3 - 
dependent terms) for solutions with (3 = 0.2, N = 1,2,3. 

Fig. 5. Energy density for iV = 3, (3 = 0.2. (A): total energy density; (B): contribution 
from (3 -independent terms; (C): GB contribution. 

Fig. 6. Yang-Mills function / for iV = 1,2,3. Solid lines: solutions with GB term 
(/3 = 0.2), dashed lines: purely EYMD solutions 

Fig. 7. Metric function W = goo (dashed lines) and exp(— 2$) (solid lines) for iV = 
1,2,3, (3 = 0.2. 

Fig. 8. Metric function a for iV = 1,2,3. Solid lines: solutions with GB term 
(P = 0.2); dashed lines: purely EYMD solutions ((3 = 0). 
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